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Abstract 

We discuss three recent developments in the calculation of Cosmic Mi- 
crowave Background (CMB) anisotropies. We begin with a discussion of the 
relativistic corrections to the Sunyaev-Zel'dovich effect. By extending the Kom- 
paneets equation to include relativistic effects, we are able to to derive simple 
analytic forms for the spectral changes due to the thermal Sunyaev-Zel'dovich 
effect. Relativistic corrections result in a small reduction in the amplitude of 
the effect over the Rayleigh-Jeans region, which for a typical cluster tempera- 
ture of 8 keV, amounts to a correction downwards to the value of the Hubble 
constant derived from combined X-ray and Sunyaev-Zel'dovich information by 
about 5 percent. Our second topic is a discussion of covariant kinetic theory 
methods for the gauge-invariant calculation of primordial CMB anisotropies. 
We present a covariant version of the Boltzmann equation, which we use to de- 
rive a set of covariant equations for the gauge-invariant variables in a Cold Dark 
Matter (CDM) model, which are independent of the background curvature and 
type of perturbation. Equations describing a particular type of perturbation 
are obtained readily from this set, as we demonstrate for scalar perturbations 
in a K = universe. By integrating the covariant Boltzmann equations along 
the line of sight, and using the instantaneous recombination approximation, we 
obtain an expression for the large scale anisotropy which agrees with an ex- 
pression derived recently by Dunsby. Finally, we discuss some recent accurate 
calculations of the CMB anisotropy in global defect models. In such models, 
the effect of vorticity generation by the causal sources proves to be significant, 
leading to a suppression of acoustic peaks. The result is that global defect 
models of structure formation may already be at variance with the growing 
volume of CMB (and large scale structure) data. 



1 Introduction 

Future satellite experiments hold the promise of obtaining the first high resolution, 
full-sky maps of the CMB. These maps will contain a wealth of cosmological informa- 
tion, allowing unprecedented accuracy in the determination of cosmological parame- 
ters, and the opportunity to distinguish between rival theories of structure formation. 
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Interpretation of the CMB data requires not only accurate calculation of the pri- 
mordial anisotropics in specific cosmological models, but also accurate modelling of 
sources of secondary anisotropies and foregrounds. In this lecture we consider three 
recent developments in the calculation of primary and secondary anisotropies. We 
begin, in Section ^, with secondary anisotropies, considering the effect of relativistic 
corrections on the spectral distortion of the CMB from Compton scattering in hot 
clusters (the Sunyaev-Zel'dovich effect). We then turn to primary anisotropies, 
presenting a summary of some work in progress on gauge-invariant calculations of 
the propagation of primordial anisotropies, using covariant kinetic theory methods, 
in Section |3]. We use a version of the covariant perturbation theory developed by Ellis 
and coworkers (see, for example, Ellis et al. @, |J). Their approach is superior to oth- 
ers, including the Bardeen formalism ||, since it deals only with physically-defined, 
gauge-invariant variables. We end with a brief review of accurate calculations by 
Pen, Seljak and Turok || |7| of the CMB power spectra in global defect theories, in 
Section [|. These calculations show that there are problems reconciling global defect 
models of structure formation with observation. 

We employ natural units c = h = G = 1, except when presenting numerical results 
in Section ||. 

2 Relativistic Corrections to the Sunyaev- 
Zel'dovich Effect 



The Sunyaev-Zel'dovich effect is concerned with the generation of CMB spectral dis- 
tortions by inverse Compton scattering in hot clusters. Non-relativistic treatments of 
this effect usually employ the Kompaneets equation |§ to determine the distortion. 
However, the Kompaneets equation does not include relativistic effects, which may be 
important for hot clusters where ksT e > lOkeV. For this reason, and because of the 
low optical depth of typical clusters, relativistic treatments of the Sunyaev-Zel'dovich 
effect usually employ a multiple scattering description of the Comptonization pro- 
cess P|-P^[|. Including relativistic effects in this procedure gives a complicated ex- 
pression for the spectral distortion, which is best handled by numerical techniques 
(see, for example, Rephaeli jl3|). 

In this section, we present an extension of the Kompaneets equation which in- 
cludes relativistic effects in a self-consistent manner. The relativistic extension to the 



Kompaneets equation has been considered independently by Stebbins JTj, [T§]. The 
extended Kompaneets equation allows the Sunyaev-Zel'dovich effect in hot clusters 
to be described analytically on the basis of a Kompaneets type equation, instead of 
the previous numerical approaches. Simple analytic forms are given for the spectral 
distortions in the limit of small optical depth, including relativistic effects to second- 
order. These are in excellent agreement with Rephaeli's numerical calculations |T3|| , 
which were based on the multiple scattering approach (truncated at one scattering). 
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This lends further support to Fabbri's observation that the Boltzmann equation can 
be applied to describe the Sunyaev-Zel'dovich effect in optically thin clusters fiCfl , 
despite claims to the contrary ||. 

The relativistic corrections to the Sunyaev-Zel'dovich effect are important in the 
calculation of the Hubble constant H by the Sunyaev-Zel'dovich route in hot clusters 
(see, for example, Lasenby and Jones |16j for a recent review of the Sunyaev-Zel'dovich 
route of determining Hq, and Saunders ]17|] for recent observations). In the Rayleigh- 
Jeans region, we find that relativistic effects lead to a small decrease in the Sunyaev- 
Zel'dovich effect, and hence a small reduction in the hitherto determined values of 



Hq, in agreement with the conclusions of Rephaeli and Yankovitch [IS 



2.1 Extending the Kompaneets Equation 

In this lecture we shall not consider effects due to the peculiar motion of the cluster 
(such effects give rise to a kinetic correction to the Sunyaev-Zel'dovich effect). For 
a comoving cluster, the CMB photon distribution function is isotropic and may be 
denoted n(uj), where u is the photon frequency. The electrons are assumed to be in 
thermal equilibrium at temperature T e , and are described by an isotropic distribution 
function f(E), where E is the electron energy. The Boltzmann equation describing 
the evolution of n(oS) may be written as fl9 l 

dn(uS) r d 3 p 



dt J (2tt) 



/ T^ d 3p ' d3k ' w r M ( x + n{uj,) ) f{E) - n{uj,) i 1 + n{uj) ) f{E,) 



(2.1; 



where W is the invariant transition amplitude for Compton scattering of a photon of 
4-momentum k^ by an electron (of charge e and mass m) with 4-momentum p^, to a 
photon momentum k' 11 and an electron momentum p ,fl [p0|] : 

= (e 2 /47r) 2 X _ 

2uu'EE' KP ' , 99 , 

1 1\ (k k /n 



X = 4m 4 - + - -V - + - - - + - . 

\K K J \K K J \K K J 

with k = ~2p fl k fl and k' = 2p fJi k' fl . In equation fl2.1p , we have assumed that electron 
degeneracy effects may be ignored. 

The electrons are described by a relativistic Fermi distribution. Since we are 
ignoring degeneracy effects, we have 

f{E) ^ e -( E -v)l k BTe_ (2.3) 
Substituting this form for f(E) into equation (|2.1|), and expanding the term in brack- 
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ets in the integrand in powers of Ax, where 



UJ 



x 



Ax = 



k B T e 
uj 1 — UJ 
k B T e 



(2.4) 
(2.5) 



gives a Fokker-Planck expansion 
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(2.6) 



(2.7) 



which does not depend on n(u). 

The calculation of the I n may be performed by expanding the integrand in powers 
ofp/m and u/m. Performing the integral over E requires an asymptotic expansion 
of the electron distribution function f(E). Evaluating the /„, we develop a (possibly 
asymptotic) expansion of dn/dt in 9 e , where 



kBT e 

m 



(2i 



In this lecture, we shall only be concerned with the lowest order relativistic cor- 
rections, and start by retaining all terms up to 0(6* e 2 ). To include all such terms 
consistently, it is necessary to retain only the first four terms in the series (2.6). A 
lengthy calculation gives the result 



dn(x) 
dt 



^NeOe^(x'j(x)), 



(2.9) 



where N e is the electron number density and the current j(x) is given by 
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(2.10) 



The zero-order term in equation (2.10) is just that term which usually appears in the 
Kompaneets equation ||. The 0(6 e ) term is the lowest-order relativistic correction 
to the current. The form of equation (|2.9| ) ensures conservation of the total number 
of photons, which is true for each order in 9 e . A similar equation has been derived in- 
dependently by Stebbins |15| using non-covariant methods, although he only includes 
terms to lowest-order in u/m which is an excellent approximation for CMB photons. 
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2.2 The Sunyaev-Zel'dovich Effect 

In this section we apply the generalised Kompaneets equation (to first-order in rel- 
ativistic corrections) to the calculation of the Sunyaev-Zel'dovich effect in optically 
thin clusters. We consider higher-order effects in the next section. 

Following the standard assumptions, we assume that the optical depth is suf- 
ficiently small that the spectral distortions are small. In this limit, we may solve 
equation (|2.9| ) iteratively. The lowest order solution is obtained by substituting the 
initial photon distribution Uq(x) into the current (eq. (2.10)). The integral over time 
is then trivial, and may be replaced by an integral along the line of sight through the 
cluster, giving 

y d ( 2 



An(x) = ^ — (x 2 j(x) , (2.11) 
x z ox K ' 

where j(x) is evaluated with n (x), and 

y = a T J N e 6 e dl, (2.12) 

where the integral is taken along the line of sight through the cluster. 

For the CMB we take the initial (undistorted) photon distribution to be Planckian 
with temperature T : 

no(x) = (2.13) 

where a = T e /T Q is the (large) ratio of electron temperature to the CMB temperature. 
Evaluating equation (|2.11|) in the limit of large a, we find the following fractional 
distortion: 

An(X) _ yXe x r„ ., ,-, „x . „/ ,„ 47 „ 42 
n(X) ~ e x - 1 



X coth(iX) - 4 + 6 e ( -10 + —X coth(iX) - —X 2 coth 2 (iX) 



— X 3 coth 3 (iX) + — (Xcoth(iX) - 3 

10 K2 J 5sinh 2 (|X)^ V2 ; 



;2.14) 



correct to first-order in relativistic effects, where 



A'^. (2.15) 

The first two terms in square brackets in equation (2.14) give the usual non- relativistic 
Sunyaev-Zel'dovich expression, while the terms proportional to 6 e are the lowest order 
relativistic correction. Equation (2.14) agrees with the result derived by Stebbins | 15fl . 
In the Rayleigh- Jeans limit (small X), we find 



n(X) 



2y(l--e e + 0(9 2 )y (2.16) 
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Figure 1: The intensity change 9 e AI/y (in units of 2(k B T ) 3 /(he) 2 ) plotted against X 
for three values of ksT e (in keV). The solid curves are calculated using the first-order 
correction to the Kompaneets equation, while the dashed lines are calculated from 
the usual Kompaneets expression. 



In Figure |l] we plot 8 e AI/y as a function of X, where 

X 3 An 
e x — 1 n 



AI = -^—, (2.17) 



is the change in spectral intensity in units of 2(/c B T ) 3 / (he) 2 . Also plotted in Figure [I] 
are the non-relativistic predictions made with the standard Kompaneets equation. 
The curves in Figure p] are for ksT e = 5, 10, and 15keV, which are the same as 
the parameters used by Rephaeli |l3j in his Figure 1. His calculations, which were 



based on the multiple scattering formalism |J and required a numerical analysis, give 
results in excellent agreement with ours, which only require the use of the simple 
expression ( |2.7| ). This suggests that there is no difficulty in principle with applying 
the Boltzmann equation to the problem of Comptonization in clusters even though 
the optical depth may be very small. Similar conclusions were reached by Fabbri [TO 



but his demonstration was restricted to low temperature clusters where relativistic 
effects are not important. 

It is clear from Figure [1] that for X<8, the relativistic corrections lead to a 
reduction in the magnitude of the intensity change, compared to the non-relativistic 
prediction. The Hubble constant is inferred from combined Sunyaev-Zel'dovich and 



X-ray data by a relation of the form fig 

# ocA/- 2 , (2.1* 
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where AI is the observed intensity change. The reduction in the magnitude of AI in 
the Rayleigh- Jeans region, for given cluster parameters, amounts to a decrease in the 
constant of proportionality in ( |2.18| ) and hence a reduction in the value of the Hubble 
constant that should be inferred. 



2.3 Higher-order Effects 

We have found that for fcsT e ^10keV the second-order relativistic effects make a 
significant contribution to the spectral distortion, while third-order effects are only 
significant for ksT e > 15 keV. 

These calculations require a straightforward extension of the method of Section [2.1| 
to include terms at 0(6^) (for second-order relativistic effects). For the calculation to 
0(# 3 ), it is necessary to retain the first six terms of the series (2.6), and to calculate 
l\ through J 6 to O(0 e 3 ). The first iteration of equation (|2.9|) for T e 3> T gives the 
following next order (in 9 e ) correction to An/n: 



2 8 y2 ' 5 

— X 3 coth 3 (iX) - -X 4 coth 4 (±X) + — _ 
5 K2 ' 5 V2 ; 30 



X 2 

+ T7T — (-2604 + 3948X coth(±X) - 1452X 2 coth 2 (|X) + 143X 3 coth 3 (±X) 

30sinlr(±X) v 



X 4 

+ jTi — rf-528 + 187Xcoth( 1 7 X) 

60sinh 4 (iX)V V2 ; 



. (2.19) 



In the Rayleigh- Jeans limit, we find 



n(X) 



2y i 1 ~ To 6 * + ^ + 0{e!) ) ■ (2 ' 20) 



In Figure we compare the spectrum of AI calculated with equation (2.14) to the 
spectrum with the correction (2.19) included, for k B T e = 5, 10 and 15keV (9 e ~ 0.01, 
0.02 and 0.03 respectively). In each case, the second-order relativistic effects are not 
significant in the Rayleigh- Jeans part of the spectrum. This is to be expected from 
inspection of equation ( j2.20| ), where the 9 2 term is clearly insignificant for the values 
of 9 e considered. For ksT e = 5keV, the second-order effects are insignificant over the 
entire spectrum. However, for ksT e ^ lOkeV, the second-order effects make a signif- 
icant contribution to the relativistic correction to the Kompaneets-based prediction 
outside the Rayleigh- Jeans region. We have verified that the third-order corrections 
are negligible over the entire spectrum for ksT e ~ lOkeV. This has been confirmed 
independently by a direct Monte-Carlo evaluation of the Boltzmann collision integral 
by Gull and Garrett |2l| . The second-order effects should be included in the analysis 



of high frequency data for hot clusters. The magnitude of the second-order correction 
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Figure 2: The intensity change 9 e AI/y (in units of 2(fcsT ) 3 /(^c) 2 ) plotted against X 
for three values of ksT e (in keV). The solid curves are calculated using the second- 
order correction to the Kompaneets equation, while the dashed lines are calculated 
from the first-order correction. 



to the Sunyaev-Zel'dovich result for the rather mild values of 6 e considered here, is 
symptomatic of the asymptotic nature of the series expansion of dn/ dt in 8 e . However, 
for the majority of clusters considered in Sunyaev-Zel'dovich analyses, the inclusion of 
the first two relativistic corrections is sufficient, particularly for experiments working 
in the Rayleigh- Jeans region of the spectrum. 



2.4 The Crossover Frequency 

The accurate determination of the crossover frequency Xo (where the thermal com- 
ponent of the spectral distortion vanishes) is essential for reliable subtraction of the 



kinematic contribution to the Sunyaev-Zel'dovich effect |13| . In Figure g we plot 
the crossover frequency as a function of fcsT e , with the first three relativistic correc- 
tions included. For fceT e < 20 keV we find that Xo is well approximated by the linear 
relation 



X ~ 3.83(1 + 1.136> e ). (2.21) 

For comparison, Rephaeli |jl3fl , found X to be approximated by X ~ 3.83(1 + 6 e ) 
in the interval k B T e = 1-50 keV, while Fabbri jTO] found X ~ 3.83(1 + 1.10 e ) for 
ksT e < 150 keV. It is clear that our calculation favours Fabbri's expression. For 
kBT e > 20keV, X calculated with the first three relativistic corrections departs from 
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k B T e /keV 



Figure 3: The crossover frequency X plotted against /csT e . The solid line is calculated 
with the inclusion of third-order corrections to the Kompaneets equation. The upper 
dotted line is a linear fit to the solid line with Xq = 3.83(1 + 1.13# e ), while the lower 
dotted line is the linear fit given by Rephaeli: X = 3.83(1 + 6 e ). 



the linear prediction ( |2.21|) . However, we do not regard this as indicative of a break- 
down of the linear approximation, since it is clear from Figure @ that the inclusion 
of higher-order terms may have a significant effect on the value of the crossover fre- 
quency. 



2.5 Conclusion 

The application of the relativistic extension of the Kompaneets equation, presented 
here, to the Sunyaev-Zel'dovich effect in optically thin clusters results in simple an- 
alytic expressions for the relativistic corrections to the usual Kompaneets based ex- 
pression for the spectral distortion. These expressions are in excellent agreement 
with the numerical calculations of Rephaeli Jl3[ for electron temperatures < lOkeV, 



providing further evidence that the low optical depth of clusters does not forbid the 
application of the Boltzmann equation to the calculation of the Sunyaev-Zel'dovich 
effect [|K|. The asymptotic nature of the series expansion of dn/dt in 9 e requires the 
inclusion of higher-order corrections to calculate the effect in hotter clusters in the 
Wien region of the spectrum. While the calculation of higher-order corrections is not 
problematic, the bad convergence properties of the series means that ultimately one 



must resort to a numerical calculation of the collision integral |^T], [22| (or employ the 
multiple scattering formalism f9j) to calculate the distortion in very hot clusters. 
Our calculations fully support the conclusions reached in Rephaeli that in- 
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eluding relativistic effects leads to a small decrease in the value of the Hubble constant, 
inferred from combined X-ray and Sunyaev-Zel'dovich information. For a cluster tem- 
perature of ~ 8 keV, the reduction in H Q due to relativistic effects is ~ 5 percent for 
measurements made in the Ray leigh- Jeans region. 

3 Covariant and Gauge- Invariant Calculations of 
CMB Anisotropics 

In this section we consider a covariant, kinetic theory approach for the gauge-invariant 
calculation of CMB anisotropies. We employ the elegant gauge-invariant perturbation 
methods of Ellis et al. || f|] (see also the recent Erice lectures by Ellis |23[ |24| and the 
contribution to this volume from Ellis and Dunsby), although we prefer to express this 
within the gauge-theoretic approach to gravity (GTG) developed recently in Cam- 
bridge, instead of conventional general relativity. GTG has been discussed in previous 



Erice lectures [25, 26], although the most detailed discussion to date is contained in 



Lasenby, Doran and Gull p7| . In GTG, gravity is introduced via gauge-fields over 
a flat (background) spacetime. The gauge-theoretic purpose of the gauge-fields is to 
ensure covariance of the theory under smooth remappings of physical events onto vec- 
tors in the background spacetime, and under local rotations of geometric quantities at 
a point in the vector space. Physical observables must be extracted from the theory in 
a gauge-invariant manner, which ensures that the background vector space does not 
play an active role in the theory. GTG differs from general relativity only in certain 
circumstances, such as the treatment of horizons P7| , solutions of general relativity 



with non-trivial topology, and the effects of quantum spin [g^, p9|] . In situations where 
none of these effects are present, such as the discussion of small departures from an 
exact Friedmann- Robert son- Walker (FRW) universe we consider here, the predictions 
of GTG coincide with general relativity. Formulating GTG in the powerful language 
of geometric algebra [|3D| allows GTG to be expressed in a coordinate-free manner, 
where all algebraic manipulations may be performed efficiently in the geometric alge- 
bra of spacetime |31[ (the STA). Our notation and conventions for GTG are the same 



as in Lasenby, Doran and Gull |27j (see also the previous Erice lectures [J25], |26|). 

The covariant and gauge-invariant treatment of cosmological perturbations pio- 
neered by Ellis et al. has many advantages over other approaches. The variables are 
covariantly-defined and relate to observable quantities in a simple manner. Further- 
more, the variables are gauge- invariant, in that they are independent of the way that 
the 'lumpy' universe is mapped onto the background FRW universe (equivalently, 
a variable is gauge-invariant in this context if it transforms covariantly under gauge 
transformations of the lumpy universe, with the gauge for the background model held 
fixed). Note that in the covariant approach, the variables are truly gauge- invariant, 
which should be contrasted to the variables in Bardeen's 'gauge-invariant' ap- 
proach, which are gauge-invariant only if the scalar, vector or tensor nature of the 
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perturbations is not altered by the remapping. A further advantage of the covariant 
approach is that the perturbation equations may be formulated without the need for 
a scalar, vector or tensor splitting of the perturbations, or the associated harmonic 
analysis, whereas these are required at the outset in the Bardeen approach. The equa- 
tions pertaining to a given type of perturbation may be obtained readily from the 
full perturbation equations by placing appropriate restrictions on the gauge-invariant 
variables. 

The covariant calculation of large scale CMB anisotropics is dealt with in the 



lecture by Ellis f32fl , where a two fluid model of the matter and radiation are used, 
and a sharp last scattering surface is assumed. On small scales, the particle-like 
nature of the radiation (and neutrinos) becomes important, as well as the compli- 
cated local physics of recombination and decoupling. An accurate treatment of small 
scale effects requires a phase-space description of the photons, and the solution of the 
Boltzmann equation describing the evolution of the distribution function. Similarly, 
a phase-space description of the neutrinos is necessary to deal properly with the heat 
flux and anisotropic stress terms in the neutrino stress-energy tensor. In this section, 
we use covariant distribution functions which we decompose covariantly into angular 
moments, as discussed by Ellis et al. |23|, |33|, |34], |35| . In these papers, the discussion 
is limited to the decomposition of the Liouville equation, appropriate for freely prop- 
agating radiation. Here we present the covariant Boltzmann equation (for Thomson 
scattering of radiation off free electrons), and its angular moments. At present, po- 
larisation is not included for simplicity. Since the generation of photon polarisation is 
sourced by anisotropy for Thompson scattering, the error resulting from the neglect 
of polarisation is small [p^j . Using this formalism, we discuss the covariant equa- 
tions describing the evolution of perturbations in a Cold Dark Matter (CDM) model. 
Our (linearised) equations are independent of scalar, vector or tensor decomposition, 
and of the spatial curvature of the background FRW model. Specialising to scalar 
perturbations in a K = universe, we obtain a set of scalar mode equations for 
gauge-invariant variables. These equations must be solved numerically to obtain the 
present day anisotropies for given initial conditions. This work will be presented else- 
where @. We end this section with a discussion of large scale anisotropies, which is 
complementary to the discussion of Ellis and Dunsby . The large scale anisotropy 
is obtained by integrating the covariant Boltzmann equation along null geodesies, and 
assuming instantaneous recombination. With the tight-coupling approximation, we 
obtain the same expression as that obtained by Ellis and Dunsby with a two-fluid 
model. 



3.1 The Gauge-Invariant Variables 

The starting point in the covariant description of cosmological perturbations is the 
choice of a velocity u. This should be chosen in a physical manner in such a way 
that if the universe is exactly FRW, the velocity reduces to that of the fundamental 
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observers. Having made a choice for u, we define a projection tensor H(a) which 
projects vectors into the space perpendicular to u: 

H(a) = a — a-uu. (3.1) 

For details of our geometric algebra notation see Lasenby et al. p5]— 1|27|| . Expressing 
the theory in geometric algebra provides frame-free and coordinate-free versions of all 
variables and operations. Using the projection ( ITT1) , we form a directional 'spatial' 
derivative, a-T* from the usual directional (covariant) derivative a-T> via 

a-VM = H[H(a)-VM], (3.2) 

where M is an arbitrary multivector. We also define a derivative a-T> which acts on 
tensors T(b,...,c) according to 

a-VT(b } ...,c) = H[H(a)-i>T(H(b), . . . , H(c))], (3.3) 

where a-VT{b, . . . , c) = a-VT(b, . . . , c) — T(a-Vb, . . . , c) — T{b, . . . , a-Vc). In an exact 
FRW universe, spatial gradients of any field M necessarily vanish by homogeneity, 
so that a ■ T>M = 0. It follows that the spatial gradient of the density T>p is a 
gauge-invariant variable. Quantities such as Dp (or, more correctly, the dimensionless 
quantity Vp/(pH), where H is a suitable average of the local Hubble parameter) will 
be small compared to unity if the universe is sufficiently close to an FRW universe. 
We shall refer to such variables as being first-order, or 0(1). In the linear theory used 
here, products of first-order quantities may be neglected in an expression, compared 
to any first-order (or zero-order) quantities. It is convenient to introduce a scale 
factor S, satisfying 

u-VS = HS, PS = 0(1), (3.4) 

so that we may remove the variation in T>p along an integral curve of u (a 'worldline') 
by defining the variable 

X^-^. (3.5) 
P 

The first-order, gauge-invariant variable X is an example of a fractional, comoving 
spatial gradient. 

The covariant derivative of u decomposes as 

a-Vu = \a-w + cr(a) + ^H(a)9 + a-uw, (3.6) 

where the vector w = u-T>u is the acceleration, the bivector w = VAu + wAu is the 
vorticity, measuring the local rotation of the worldlines, and the scalar 9 = T>-u = 3H 
measures the volume expansion rate. The traceless symmetric tensor a (a) measures 
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the shear of the worldlines, and is orthogonal to u: cr(u) = 0. We shall refer to 
tensors orthogonal to u as being spatial tensors. The vorticity, acceleration and 
shear are all first-order gauge- invariant variables (they vanish in an FRW universe). 
However, like p, the volume expansion rate 9 is zero-order, but we obtain a first-order 
gauge-invariant variable Z by taking the spatial gradient, Z = ST>9. 
The matter stress-energy tensor T{a) decomposes as 

T(a) = pa-uu + qa-u + uq-a — pH(a) + 7r(a), (3.7) 

where the scalar p = T(u) ■ u is the density of matter (measured by a comoving 
observer), the vector q = H[T(u)} is the energy (or heat) flux, the scalar p = —H(d a )- 
T(a)/3 is the isotropic pressure, and the traceless, symmetric, spatial tensor ir(a) = 
H[TH(a)] + pH(a) is the anisotropic stress. In an exact FRW universe, the stress- 
energy tensor is forced to take the ideal fluid form, hence q and 7c(a) are first-order 
variables. A first-order, gauge-invariant variable may be derived from the pressure 
by taking a spatial gradient. 

The final gauge-invariant variables we shall make use of derive from the Weyl ten- 
sor W(B), which is a bivector- valued linear function of a bivector B. It is convenient 
to split the Weyl tensor into electric and magnetic parts, in analogy with the split of 
the Faraday bivector in electromagnetism. We define the electric part by 

S(a) = -w-W(aAw), (3.8) 

and the magnetic part by 

B{a) = -iuAW(aAu), (3.9) 

where i is the unit pseudoscalar for spacetime (which performs the duality opera- 
tion in the STA). Both S(a) and B(a) are traceless, symmetric spatial tensors. The 
definitions ( |3.8| ) and ( |3.9| ) combine together to give the neat relation 

W(aAu)u = 8(a) + iB(a). (3.10) 

Since the Weyl tensor vanishes in an FRW universe, the electric and magnetic parts 
of the Weyl tensor are first-order gauge-invariant variables. The equations of motion 
for the gauge-invariant variables may be derived from the Ricci identity, the Bianchi 
identity, and the contracted Bianchi identity. These equations are the GTG equiva- 



lents of the equations given by Ellis |32| in his lecture, and will be given in linearised 
from in Section |3.4| . The equations split into propagation equations (those involv- 
ing the derivative along a worldline u-V), and constraint equations which involve 
only spatial derivatives. The constraint equations restrict the specification of initial 
data on a hypersurface (the data must satisfy the constraints), while the propagation 
equations allow one to propagate the data off the surface. 
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3.2 The Boltzmann Equation 

The photon distribution function is a scalar valued function of position x and covari- 
ant momentum p, denoted by f^(x,p). Physically, an observer sees f 1 (x,p)dV\d 3 p\ 
photons in a volume dV and a 3-momentum volume \d 3 p\. It is convenient to split 
the photon momentum p as 

p = E(U + e), (3.11) 

where E = puis the photon energy and e is a unit spacelike vector e 2 = — 1 orthogonal 
to u. The photon distribution function may then be written as / 7 (x, E,e). We will 
usually leave the x dependence implicit, writing f 1 {E,e). The stress-energy tensor 
for the photons may be written as 

T 7 (a) = J dEdttEf^(E,e)pp-a, (3.12) 

where dfl denotes an integral over solid angles. 

Following Ellis et al. |35j, we describe the angular dependence of f^(E,e) by 
expanding in covariant, scalar- valued tensors if (ai, a 2 , ■ ■ ■ , a{): 



oo 

/ 7 (£,e)=£if (e,e,...,e). (3.13) 



1=0 



The tensors if (ax, a 2 , ... , ai) have an (implicit) dependence on x and E and satisfy 
the properties: 

if (ax,... 

, Ctj , . . . , flfc , . . . , d[ ) = if K... 
F«(«,a 2 ,...,a,) = (3.14) 
F«(^,6,a 3 ,...,o,) = 0. 

The action of the Liouville operator C on if (e, . . . , e), evaluates to 

£if (e, . . . , e) = ^F«(e, . . . , e)d x E + p-VF®(e, . . . , e) + lF®(p-Ve, e, . . . , e), 

(3.15) 

where A is the affine parameter to the photon path (p = hr l (d\x) with h(a) the GTG 
position-gauge field). It follows that the action of the Liouville operator on f^(E,e) 
may be written in the form 

oo 

Cf,(E, e) = X> B if (e, . . . , e)d\E + p-VF«\e, . . . , e) + lF^(p-Ve, e, . . . , e)]. 

'=o (3.16) 

For I > the tensors if (ax, . . . , a/) are first-order quantities. However, to zero-order, 
we have 

p-Ve = \6Eu, (3.17) 
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so that the final term in ( |3.16| ) only contributes at second-order, and may be ignored 
in the linearised calculation considered here. 

Since neutrinos are essentially collisionless over the period of interest, their dis- 
tribution function f u (x,p) evolves according to the Liouville equation Cf v (x,p) = 0. 
The photons, however, are in close contact with the baryons through Thomson scat- 
tering off free electrons (which are at rest in the frame of the baryons because of 
Coulomb interactions). In the presence of collisions, the photon distribution function 
evolves according to the Boltzmann equation 

Cf J (x,p) = C, (3.18) 

where the collision term approximates to 

C = n e a T p-u b (f + (x,p) - f(x,p)), (3.19) 

with Ub denoting the baryon velocity, ot the Thomson cross section, and f+(x,p) 
describes scattering into the phase space element under consideration. Evaluating 
f+(x,p) to first-order, multiplying by E 2 and integrating over energy, we find |2j 

J dE E 2 Cf 1 (E 1 e) = YQ^ aTTle (3 ~ 4e '^)P7 + e " TC -y( e )J 

POD 

-a T n e dEE 3 UE,e), (3.20) 
Jo 

where v b = u b — u is the first-order relative velocity of the baryons. 

Equation (3.20) may itself be expanded in symmetric, trace-free tensors orthogonal 
to u (like Fjp(ai, . . . ,ai)). The / = 0, 1 and 2 terms require a little care in the 
linearised calculation since F^ is a zero-order quantity. The result is @: 
for I = 

u-V Pl + f#p 7 +V-q 1 = 0, (3.21) 

for / = 1 

u-Vq 1 + |#g 7 + 7r 7 (P) + |u>p 7 - \Vp^ = a T n e (|p 7 ^ - <? 7 ) , (3.22) 

for I = 2 

ar(u-P7r 7 (a 2 )) + §0ai-7r 7 (a 2 ) + J^(V,a 1 ,a 2 ) - l(a 1 -(a 2 -Vq J ) + a 2 -{a 1 -Vq 1 ) 

- ^H(a 1 )-H(a 2 )V-q 1 J - j^a 1 -a(a 2 )p J = -j^a T n e ai-n 1 (a 2 ), (3.23) 

and for I > 3 

u-Vj®(a u . . . , a{) + pJ^iax, . . . , a,) + J^\f>, a 1; . . . , a,) 
~~ {2i+\)i\ ( a i-'DJy~ 1 \a 2 , ...,ai)- ■^^J^ l ~ 1 \V,a 1 , . . . , ai- 2 )H{ai-i)-H{ai) + perms J 

= -a T n e jW(a 1 ,...,a l ). (3.24) 
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The scalar- valued, traceless, symmetric, spatial tensor jW(ai, 

4tt(-2) / (/!) 2 f°° 



ai) is denned by 



"i v«i,---,^- ( 2 / + i)(2/)!y 
With this definition, we have 



P 7 = ^ 0) 



dEE 3 F®{ ai ,...,ai). 



q 1 = d a j\ l \a) 
7r 7 (a) = d b A 2) (b,a). 



(3.25) 



(3.26) 
(3.27) 
(3.28) 



It is not hard to show that the combination 



,ai- 2 )H(ai-i)-H(ai) + perms 



7 

(3.29) 



is a trace-free, symmetric tensor, orthogonal to u, as required. Note that the linearised 
equations show a coupling between the I — 1, I and Z + 1 terms, whereas the exact 
equations also show coupling between the I — 2 and I + 2 terms. The terms involving 
J~~ 2 ^{0'i, • • • , 0*1-2) an d J^ +2 )(ai, . . . , cfy+2) also involve the shear tensor a (a). Ellis 
has remarked that the exact result that if the hierarchy of Jjp tensors truncates 
after a finite number of terms then the shear must vanish, is missed in the linearised 
approach. This is certainly true, however it is not problematic for the calculation 
of CMB anisotropies, since it is never claimed that the expansion of the photon 
distribution function truncates. Instead, the series is truncated (with suitable care 
to avoid reflection of power) after a finite number of terms for numerical convenience 
(see, for example, Ma and Bertshinger |37]). The truncation is chosen to be high 
enough up the series that it has negligible effect on the Jjp for the range of I of 
interest. 



3.3 Baryons and Dark Matter 

Over the epoch of interest here, the baryons and electrons are non-relativistic, and we 
assume that they may be described as an ideal fluid, with energy density pb (measured 
in the rest frame of the baryons), pressure p&, and covariant- velocity u + Vb, where Vb 
is a first-order quantity. The linearised baryonic stress-energy tensor evaluates to 

%{a) = pbO-uu — p b H(a) + (p b +pb)(a-uv b + a-v b u), (3.30) 

which shows that there is a heat-flux (p b + PbjVb measured by an observer moving 
with velocity u. To find the equations of motion for pb and v b , we make use of the 
fact that the baryons (including electrons) and photons interact non-gravitationally 
only with themselves, so that 

t / (V) + f b (V) = 0. (3.31) 
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Using the moment equations of the previous section, we find a propagation equation 
for p b : 

u-Vp b +(p b + p b )9 + (p b +p b )V-v b = 0, (3.32) 

and a propagation equation for v b \ 

(Pb +p b )(u-Vv b + w) + \{p b + p b )8v b + u-Vp b - Vp b + a T n e (fp 7 f& - ? 7 ) = 0, 

(3.33) 

which must be supplemented by an equation of state linking p b and p b . Note that 
equation ( |3.32| ) implies that there is no energy exchange between the radiation and 
the baryon plasma in the frame in which the baryon plasma is at rest, while equa- 
tion ( |3.33j ) shows that there is momentum exchange due to dipole anisotropy of the 
CMB in the frame of the plasma. This behaviour is to be expected since we have 
assumed that the baryon-photon interaction approximates to Thomson scattering in 
the baryon rest frame, for which there will be no energy transfer, but there will be 
momentum transfer if the photon distribution function is not isotropic. 

We will only consider cold dark matter (CDM) here, which may be described as a 
pressureless ideal fluid. Hot dark matter (HDM) would require a distribution function 
description which for massive particles greatly increases the numerical complexity of 
the problem. (Both CDM and HDM are considered, for example, in the gauge- 
dependent treatment in Ma and Bertshinger |j37j|). The CDM has energy density p c 
in its rest frame, which has covariant velocity u + v c , with v c a first-order quantity. 
The CDM interacts with other species through gravity alone, so the equations of 
motion are 

u-Vp c + p c 9 + p c V-v c = (3.34) 
u-Vv c + l6v c + w = 0. (3.35) 

It is very convenient to use the CDM velocity to define the fundamental velocity u 
(v c = 0). Since the CDM is pressureless, it moves geodesically so that the acceleration 
w will vanish for this choice of u. The discussion in the next two sections assumes 
that this choice has been made. 

3.4 General Equations for the Evolution of CMB Anisotropies 

As Ellis and Bruni || have shown, the natural variables with which to describe 
cosmological perturbations are the comoving fractional spatial gradients of the energy 
densities, where i labels the particle species. The equations of motion for these 
variables follow from taking the spatial gradients of the evolution equations for the 
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pi. This procedure yields the linearised equations: 

u-VX u = --Z - —VV-q v (3.36) 

3 Pu 

4 S - - 

u-VX~ = — Z VV-q, ( (3.37) 

3 p 7 

n-D^ c = -Z (3.38) 

u-VX b = -(l + m)(Z+ SVV-v b ) - S6^ + ^ 6 ^. (3.39) 
V PbJ Pb pb 

An equation of motion for Z may be obtained by taking the spatial gradient of the 
propagation equation for 6 (the Raychaudhuri equation, which follows from the Ricci 
identity) to obtain 

u-VZ = -\QZ - k \p 1 X 1 + p v X v + \p c X c + \p b X b + \SVp b 



(3.40) 

Note that these equations do not close due to the presence of v b , g 7 and q v . These 
equations, along with the kinetic theory equations given in the previous section and 
the equations for the gauge-invariant variables, discussed in Section |3~T| , form a com- 
plete description of the evolution of CMB anisotropies. For convenience, we group 
these equations together in this section. 

The baryon peculiar velocity v b evolves according to 

(Pb + p b )u-Vv b + \{p b + Pb)9v b + v h u-Vpb - Vpb + cr T n e (§p 7 t> 6 - g 7 ) = 0. 

(3.41) 

For the neutrinos, we have the moment equations 

U'Vq v + ±0q v + Tt v {p)-\Vp v = Q, (3.42) 
a 1 -(u-7>k 1/ (a 2 )) + %6ai-ir u (a 2 ) + JjP(f>, a%, a 2 ) - \{ai-(a 2 -Vq v ) + a 2 -(ai-T>q„) 

- lH{a x )-H{a 2 )V-q v ) - ^ ai -a(a 2 )p u = 0, (3.43) 

and for I > 3 

u ■ VJl l) (a x , . . . , at) + §0 4 l) (a x , . . . , a t ) + [V, a x , . . . , ai) 

- wtw( avf>J »~ l) ( a2 > ■ ■ ■ ' ai) ~ §TT) J t 1) ( f) ^ °i' • • • > ai-zMaj-i) •#(<!,) + perms) 

= 0. (3.44) 

For the photons, we have 

u-Vq 1 + |#g 7 + TTryiV) - §£>p 7 = cr T n e (fp 7 w ft ~ Q-y) , (3.45) 
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ai-(u-t>7c J (a 2 )) + |6>a r 7r 7 (a 2 ) + Jff\T>, a 1; a 2 ) - ^a r (a 2 -'Dg 7 ) + a 2 -(a r Z>g 7 ) 

- |if(ai)-if(a 2 )P-g 7 ) - ^ar<7(a 2 )p 7 = -^<T T n e ar7r 7 (a 2 ), (3.46) 

and for I > 3 

u-X> j«(ai, . . . , a,) + tflJ^K . . . , a,) + J 7 ' +1) (P, «i, • • • , aj) 
~ (2i+i)ii (arT>J^- 1 \a 2 , . . . , a,) - ji 1 ' 1 )^, a u . . . , a l ^ 2 )H(a l _ l )-H(ai) + perms) 

= -<r T n e J®(a 1 ,...,ai). (3.47) 

The remaining equations are the propagation and constraint equations for the 
kinematic variables and the electric and magnetic parts of the Weyl tensor. The 
constraint equations are 

B(a) - \[(iaAuAV)-m + iuAVA(a-m)\ + \[iuAVAo(a) + o(iuAaAV)\ = 

(3.48) 

B(D) = -\k[(p + p)iuAw + iuAVAq] (3.49) 
£(V) = \k[2V P + 29q + 3tt(P)] (3.50) 
\V-w-a(V) + \V6 = -nq (3.51) 
V-(iuAw) = 0, (3.52) 

where p and p are the total energy density and pressure measured by an observer 
moving with velocity u: 

P = Pv + P 7 + Pb + Pc (3.53) 
p = \p v + |p 7 + p b , (3.54) 

and q and re {a) are the total heat flux and anisotropic stress in the frame defined by 
u: 

q = q u + q y + (p b + p b )v b (3.55) 
7r(a) = 7r u (a) + 7r 7 (a). (3.56) 

Note that p and p are independent of the choice of u to first-order, and that there is 
no CDM contribution to q since we have chosen u equal to the CDM velocity. 
The propagation equations are 

- u-V£(a) - 68(a) - l B (a) = ^re[6(p + p)cr(a) + 3(a-f>q + V(a-q)) 

- 6u-T>Tc(a) - 29n(a) - 2H(a)V-q] (3.57) 

- u-VB(a) - 9B(a) + l £ (a) = -\n[iu AV An (a) + n(iuAaAV)] 

(3.58) 
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u-Va(a) + \9a{a) + £(a) + \mx(a) = (3.59) 

u-Vm+l6m = 0, (3.60) 
where Xg(a) is a trace-free, symmetric, spatial tensor defined by 

J B (a) = §[mAXM£(a) + 6(mAoAP)], (3.61) 

where 

B(iuAaAV) = H[H(b)-VB{iuAaAd b )}, (3.62) 

and 2jp(a) is defined similarly. There is some redundancy in these equations. For 
example, equations ( |3.48| ), ( |3.51| ) and the integrability condition 

VAV9 = -zuu-Ve, (3.63) 



imply the constraint Q3.49Q . Similarly, equation ( 3.48 ) along with the propagation 



equations ( 13.591 ) and (|3.60|) imply the propagation equation ( |3.58j) . It follows that 



B(a) may be eliminated from the full set of equations by making use of the con- 
straint (|3.48p . This turns out to be a necessary step when the covariant first-order 



quantities are harmonically expanded for scalar, vector or tensor perturbations. 

Some comments are in order regarding the equations presented in this section. The 
equations are both covariant and gauge- invariant (in the sense of being independent of 
any map between the lumpy universe and the FRW model). Gauge-invariance ensures 
that the gauge-problems that have plagued other approaches are not present, while 
covariance ensures that we are working with quantities which are straightforward to 
interpret physically. The equations describe scalar, vector and tensor perturbations 
in a unified manner, and are independent of any harmonic analysis into spatial modes. 
Furthermore, we have not had to specify the background FRW model yet (there is 
an implicit assumption that the universe is approximately FRW when the first-order, 
covariant variables are constructed). 



3.5 Scalar Perturbations in a K = Universe 

In this section we reduce the general equations of the previous section to a set of 
equations for scalar-valued, first-order gauge-invariant variables describing the evolu- 
tion of the density inhomogeneities and the CMB anisotropy for scalar perturbations 
in a universe which is approximately a K = FRW model. These equations thus 
describe the standard CDM model in a covariant and gauge-invariant manner. The 
equations split into a set of algebraic relations (from the constraint equations) and 
propagation equations for scalar variables. The moment equations for the photon and 
neutrino distributions (for I > 3) are equivalent to those given elsewhere |37H , where 
a Fourier expansion of the spatial variation, and a Legendre expansion of the angular 
dependence of the CMB anisotropy are made. 
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Scalar perturbations may be characterised in a covariant manner by demanding 
that the magnetic part of the Weyl tensor and the vorticity both vanish identically. 
The vanishing of the vorticity ensures that u is a hypersurface orthogonal vector 
field. The perturbation equations place strong restrictions on the remaining non- 
zero variables, which may be satisfied by constructing the vector and tensor variables 
from spatial derivatives of the scalar (harmonic) eigenf unctions Qk of the generalised 
Helmholtz equation: 



V 2 Q k = %>Q k , (3.64) 



which are constructed to satisfy 



u-VQ k = 0(1). (3.65) 
From the Q k we form a scalar valued tensor Q^\a): 

Q£\a) = la-VQ k , (3.66) 

which has the properties 

Q k 1 \u) = 0, u-VQ k 1) (a) = 0(1). (3.67) 

We then define scalar- valued tensors Q k \ai, . . . , a{) by the recursion formula (for 
1>1) 

• • • , a t ) = if (arPQ^K • • • , a,) - ±±H( ai ) ■ H (a^Q^ (V , a 3 , . . . , a,) 

+ perms). (3.68) 

These tensors satisfy the properties 

) = Qt'(a 1 ,... 

i *^fc j • • • > j • • • j O'l) 

Q k l \u,a 2 , • • • ,a t ) = 
Q k l \d b ,b,a 3 , ...,ai) = 

u.PQg ) (ai,...,o,) = 0(l), (3.69) 

which are readily proved by induction. We shall also make use of the vector Q k = 
d a Qk\ a ) an d the vector-valued tensor Qk{a) = dbQ k 2 \a,b). 

The definitions above are independent of the spatial curvature K of the back- 
ground model. However, the differential properties of the Q k (ai, . . . , Oj) are depen- 
dent on the value of K. Some of these properties are listed in the appendix to Bruni et 
al. |38] The only new result we require is 

Qf (V, a 2 , . . . , ai ) = ^ifQ?^ (aa, . . . , a,), (3.70) 
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which is valid for K = only. 

We separate out the spatial and temporal dependence of the covariant, first-order 
variables by expanding in the harmonics derived from the Q k in the following manner: 



k 


Z 


= E k iz k Qt ] 

k 


(3.71) 


S(a) = ^2(l) 2 ^ k Q k (a), 

k 


a (a) 


= E (!) °kQk(a) 

k 


(3.72) 


qu = pu^qukQ[ V \ 

k 


<h 


= Pi E 1-ykQk 

k 


(3.73) 


n u {a) = p u ^Z^ukQki.a), 
k 


7r 7 (a) 


- p 7 E 7r 7fc ( 3fc( a )- 

k 


(3.74) 



We also expand v b as 



v b = J2v k Q^. (3.75) 



Finally, we assume that the higher-order moments of the (energy-integrated) neutrino 
and photon distribution functions may also be expanded in harmonics. By consid- 
ering the zero-order form of the scalar harmonics Q k , it is straightforward to show 
that this assumption is equivalent to that which is usually made |37j (the Fourier 
components of the distribution function are axisymmetric about the wavevector k). 
The procedure we have adopted here, of expanding the angular dependence of the dis- 
tribution function in symmetric, traceless, covariant spatial tensors, which are then 
expanded in the appropriate harmonic tensors, derived from the scalar harmonics, 
should be compared to the usual approach of Fourier expanding (for K = 0) all vari- 
ables, and then performing a Legendre expansion of the angular dependence of the 
distribution function modes about the wavevector k. Although ultimately equivalent, 
the approach taken here has the advantage that it goes over unchanged to the case 
of tensor perturbations (although it is now Q k 2 '(ai, a 2 ) which satisfies the generalised 
Helmholtz equation). For later convenience, we expand the angular moments of the 
distribution function in the form 



J«\a u . . . , a,) = p 7 E J?kQk(au ...,a t ). (3.77) 



4 l \a u ...,a l ) = p u J2 4M>i, • • • > at) (3-76) 

k 

K0r>(0 ( 



It is now a simple matter to substitute the harmonic expansions of the covariant 
variables into the general equations of Section (with B(a) and w set to zero), to 
obtain equations for the scalar expansion coefficients. We assume that the variations 
in baryon pressure pi, due to entropy variations are insignificant compared to those 
due to variations in p b . Accordingly, we write 

V Pb = c 2 s V Pb , (3.78) 
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where c s is the adiabatic sound speed in the baryon/electron fluid (this is different 
from the sound speed in the coupled baryon/photon plasma). 

With this assumption we obtain the following equations for scalar perturbations 
in a K = universe: the spatial gradients of the densities give 

Xvk — — |f %k — f <2Vfc (3.79) 

X-yk = — |f Zk — f <?7fc (3.80) 

X bk = - (l + a) |(z fe + Vk ) + c 2) 9X bk (3.81) 
^cfc = — f %ki (3.82) 
(an overdot on a scalar denotes u-V), the spatial gradient of 6 gives 



k 

s' 



\Z k + |f dZ k + i/t(2p 7 A' 7fc + 2p v X vk + (1 + Sc^p,^ + p c ^ cfc ) = 0, 
the heat fluxes give 



(3.83) 



Qvk + §5 7r ^fe ~~ |f <^fc — (3.84) 

Qjk + ff^Tfe — if^yfc = ^T^e (|^fc — <?7A;) (3.85) 

(l + J) (t* + 1(1 - 3c^ fc ) - c%X hk = - P f b O T n e (lv k - q lk ) , 

(3.86) 

the remaining moment equations are, for / > 3 

• , 3 k t(3) 2fc^8fe^. n /o o7\ 

^ 7 fe + |f J$ - if ?7fc _ ilf = -B°"Tn e 7r 7A: (3.88) 

1(0 , fe + tC + 1) l_ j(J 1)\ _ A /q OQ\ 

^i/fc 5 \2l+i J vk 2l+i J vk ) — u l^o.oy; 
J jk 5 V2«+l J 7fc 22+1 7^ ) — "T^e-J^k, [0.i)V) 

the electric part of the Weyl tensor gives 

- (f ) 2 + |0$ fc ) = 5«| ((p + p)<T fc + p^^jfc + p 7 g 7fc + (p 6 + Pb)v k ) 

+ \K6(p v -K vk + p 7 7T 7A: ) - \K(p v -K vk + p 7 7T 7fc ), (3.91) 

and the shear tensor gives 

f (<r fc + + (f ) 2 + ^(P^^fe + P 7 vr 7fc ) = 0. (3.92) 

Finally, the harmonic expansions of the constraint equations give 

2 (f ) ~~ K^(p v (X vk + lT uk ) + p 7 (A' 7 fc + 7T 7fe ) + p c X ck + PftAfftfc) 

- K0(p v q v h + P-yQ-yk + (pb + Pb)Vk) = 0, (3.93) 
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and 

I (I) ( Zk ~ ak ) + K (p^k + Pyqyk + ifib + Pb)vk) = 0. (3.94) 

It is straightforward to show that these constraint equations are preserved by the 
propagation equations. 

These gauge-invariant equations should be compared to the (gauge-dependent) 
equations used in most calculations of the CMB anisotropy (see, for example Ma and 
Bertshinger |37|). Note that the equations for J^l and for / > 3 are equivalent to 
those usually found in the literature since the moments of the distribution function 
in such gauge-dependent calculations are actually gauge- invariant for / > 3. 



3.6 The CMB temperature anisotropy 

The equations of the previous section may be solved numerically for appropriate initial 
conditions [[| to give the J® at the present epoch. These moments fully describe the 
CMB temperature anisotropy. Denote the full sky average of the CMB temperature 
by T , and the (gauge- invariant) temperature difference from the mean along a spatial 
direction e by 5T(e). Then we have [|23| 

(3.95) 



(T + 5T{e)Y EZo I dE £ 3 T« (e, . . . , e) 



so that, to first-order, 



T$ J dE E 3 F^ 



^ = -Y,fdEE*F«\e,...,e). (3.96) 



Pi 1=1 

Recalling the definition ( |3.25| ), we find 

5T(e) _ 3q r e 15 7r 7 (e)-e 1 » (2l + l)(2Q! (l) 

T 4 p 7 + 8 p 7 4p 7 ^ (-2)'(J!)* f 1 ' 



(3.97) 



The right-hand side of ( |3.97p is the covariant harmonic expansion of the temperature 
anisotropy. The temperature anisotropy may be expanded in the more familiar spher- 
ical harmonics by introducing a covariant triad, orthogonal to u, at the observation 
point. We now have all the ingredients needed to calculate the CMB anisotropy for 
scalar perturbations in a flat CDM model, for given initial conditions. The numerical 
calculation of the anisotropy for such a model are underway, and the results will be 
reported elsewhere ||. Our results should confirm the calculations of other groups, 
once the gauge-invariant information is extracted from their results (the CMB power 
spectrum, Ci, for / > 0). 
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3.7 Large Scale Anisotropies 

On large angular scales, the dominant contributions to the anisotropy may be ex- 
tracted analytically without recourse to a full numerical treatment. In this section, 
we show how this can be done within the framework developed above. Our gauge- 
invariant treatment of large-scale anisotropies is complementary to the discussion by 
Ellis and Dunsby [E2] elsewhere in this volume, who use a two-fluid model for the 



matter and radiation after recombination, rather than the kinetic theory approach 
adopted here. 

For semi-analytic work, it is convenient to use the Boltzmann equation in the form 
of equation (3.20). Defining 

S T (e) ee (3.98) 
-to 

where ST(e) is the gauge-invariant temperature fluctuation from the mean (see equa- 
tion (|3.95| )), we can write (3.20) in the form 

1 p' 
6 T {e)' + a T n e S T (e) = e-a(e) + e-w - -0(1 + 45 T (e)) - f^(l + 4<5 T (e)) 

3 4p 7 

(3.99) 

where a prime denotes differentiation with respect to the parameter A along the 
geodesic, where, for this section, we take (u + e)-T>\ = 1. Equation (3.99) is correct 
to first-order. Note that in writing (3.99) we have not made a physical choice for the 
velocity u yet, so we allow the possibility that w does not vanish. (In the previous sec- 
tion we took u to coincide with the velocity of the CDM, but we relax that restriction 
here.) Writing the Boltzmann equation in the form (3.99) is useful since the equa- 
tion may be integrated along the geodesic to determine the temperature anisotropy 
along the given direction on the sky. Before doing this, it is convenient to eliminate 
9 from (3.99). In this section we shall consider a model with interacting baryons and 
radiation only, so that we describe the same situation as Ellis and Dunsby elsewhere 




in this volume JSfl (see also Dunsby P9|). In this case, it is convenient to define u to 
be equal to the velocity of the baryons, so that Vb = 0. Using ( |3.32j ), and neglecting 
baryon pressure, we find 

6 = - P -± + ^ (3.100) 
Pb Pb 



and using ( |3.33|) , we may write the acceleration as 



w = <TTn e —, (3.101) 
Pb 
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which shows that after recombination, the acceleration is negligible, particularly in a 
universe which is matter dominated at recombination. Substituting into (3.99), we 
find 



S T (e)' + a T n e S T (e) 



e-a(e) — 



1 e-Vp b 



Pb 



P'b P 7 , „ „ ( e ^7 , 3 e-7r 7 (e) 

h a T n e 1 

3p b 4p 7 V p b 



16 



Pi 



(3.102) 



This equation may be integrated up the null geodesic from some point in the distant 
past (where A = A») to the reception point R (where A = Xr). Introducing the optical 
depth k(X) along the line of sight, defined by 



k(X) = 

we find that ( |3.102| ) integrates to give 



n e <jT dX, 



(3.103) 



(*r(e)) 



R 



e-q-y 3 e-ir~(e) 1 , 1 , 

— + T7i — - - hi p b + - In p 7 

Pb 16 p 7 3 4 



+ e~ 



e-a(e) 



l e-Vp b 
3 pb 



dX, (3.104) 



where we have integrated by parts, we have assumed that terms evaluated at Aj 
are negligible due to the large optical depth, and we have neglected a direction- 
independent (monopole) term which must be cancelled by other terms in the integral. 
The notation (A)r denotes the quantity A evaluated at the point R. An expression 
similar to (3.104) forms the starting point of the 'line of sight integration approach' to 
calculating CMB anisotropies [[4(J, but note that (3.104) is true for scalar, vector and 
tensor modes in almost FRW universes with any value of K, and has not assumed any 
Fourier decomposition. Furthermore, the expression only involves physically-defined 
gauge-invariant variables, such as the fractional temperature fluctuation from the 
mean &r(e). 

The quantity —n'e~ K defines the visibility function, which is the probability den- 
sity that a photon was last scattered at A. The visibility function peaks at a red- 
shift z ~ 1100 and has a dispersion ~ 70 in redshift |fjl| . It follows that on an- 
gular scales larger than 8'f2 _1//2 , the visibility function may be approximated by a 
5-function, whose support defines the last scattering surface. With this approxima- 
tion, e~ K = before last scattering (tight-coupling), and e~ K = 1 after last scattering 
(free-streaming), so that equation (3.104) reduces to 



(*r(e)) 



e-q 1 3 e-7r 7 (e) 
Pb 16 p 7 



-lnp 6 + -hip^ + 



e-er(e) — 



1 e-Vpb 



Pb 



dX, 



(3.105) 
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where A is the point where the geodesic intersects the last scattering surface. Taking 
the difference between two directions on the sky, we find the following expression for 
the gauge-invariant temperature difference AT 



/AT 



V T 



/ R 



l An v l An v A(e-g 7 ) E 3 A(e-7r 7 (e)) E 
-A(lnp 7 ) E - -A(lnp 6 ) E + '— + 



4 



Pb 
+ A 



16 



Pi 



le-Vp b 

e-a(e) - - 

3 Pb 



where, for example, 



A(lnp 7 )£ = (lnp 7 ) A - (lnp 7 ) B , 



dX, (3.106) 



(3.107) 



with B the point of intersection of the other null geodesic with the last scattering 
surface, and 



A 



dX 



( 



dX- 



Afl 

( 

Ab 



dX. 



(3.108) 



The first term on the right-hand side of (3.106) is just the fractional difference in 
radiation temperature between the points of emission A and B on the last scattering 
surface, while the second term describes the effects of inhomogeneity in the baryon 
density. It is convenient to write the sum of these two terms as a line integral in 
the last scattering surface. Denoting by dx cov the covariant element of length for an 
arbitrary path connecting A and B, we have 



1 1 

-A(lnp 7 ) B - -A(lnp b ) E 



f 


T>p 7 


IB 


4p 7 


f 


Yr)p 7 


IB 





Vpb 
3pb 

_ 

3p6 



+ u 



'u-Vp 7 u-Vpb 



4pn 



3p& 



B 



[S^b ~ \uD-v^\ -dx cov , 



where 



S. 



■yb 



T>p 7 Vp b 



v 1 = 



3g 7 



• dx cov 
(3.109) 

(3.110) 



4p 7 3p b 4p 7 

The covariant vector 5 7 b is a covariant measure of the entropy perturbations (see, for 
example, Bruni et al. [3~8"]), and t> 7 is the effective radiation velocity (relative to the 
baryon frame). If we take the integration path to lie entirely in the last scattering 
surface, u-dx cov = 0(1), so that (3.109) reduces to 



1 1 [A 

-A(lnp 7 ) jB - -A(\np b ) E = J b S~ fb -dx c 



(3.111) 
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which shows that the contribution of the first two terms in (3.106) depends only on the 
entropy perturbations on the last scattering surface. In particular, for perturbations 
which are exactly adiabatic at last scattering, there is no contribution from these 
terms. The third term in (3.106) is a Doppler term which is negligible for a universe 
which is matter dominated at recombination, and the fourth arises from viscous effects 
in the photons. Before recombination, both g 7 and 7r 7 are first-order in the photon 
mean free time {n e aT)~ l (recall that we have defined u to coincide with the baryon 
velocity), so that in the tight-coupling/instantaneous recombination approximation 
these terms may be ignored. It follows that for perturbations which are exactly 
adiabatic at last scattering, only the integral term on the right-hand side of (3.106) 
remains. 

An expression similar to (3.106) is derived elsewhere in this volume by Ellis and 
Dunsby [|3~2]] , although they used a two fluid approach rather than the kinetic theory 
approach adopted here. To make contact with their work, we consider a universe 
which is already matter dominated at recombination. In this case, the baryon velocity 
coincides with the velocity of the energy frame (the frame in which the heat flux 
vanishes) which was used by Ellis and Dunsby, and we may replace pb by the total 
density p in the integral in ( |3.105|) . Furthermore, we ignore the terms involving g 7 



and 7Try , CIS discussed above, to obtain 



11/" 
( s T(e)) R = -(lnp 7 )A - -(hipbU + J x 



A.i 



e-cr e 



le-Vp 
3~p~ 



dX. 



(3.112) 



For ease of comparison we also give the corresponding expression from Ellis and 
Dunsby [^| in our notation: 



(*r(e) 







)n = A + 


[ Xr 




J\ A 



e-a(e) 



le-Vp 
3~p~ 



dX, 



(3.113) 



where 



A 



e-(S. 



■yb 



^V-v^) dX. 



(3.114) 



By making use of the rearrangement in (3.109), and noting that dx cov = (e + u)dX 
along the geodesic, we find that A reduces to the first two terms on the right-hand side 
of (|3.112|) (after subtraction of a monopole term). It follows that the expression for the 
temperature anisotropy in Ellis and Dunsby is in agreement with that derived here, 
under the assumptions outlined above. Note that the part of A which is significant 
observationally is determined only by the entropy perturbation at last scattering, and 
so this term will vanish for perturbations which are adiabatic at last scattering, even 
though the perturbations do not remain adiabatic after decoupling. 
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3.8 Conclusion 



The gauge-invariant calculation of CMB anisotropics may be performed in a fully co- 
variant manner with the methods outlined in this section. The equations of Section p.4| 
are a complete description of the evolution of perturbations in a CDM universe. They 
include all types of perturbation (scalar, vector and tensor) implicitly and are inde- 
pendent of the curvature of the background FRW model. We have demonstrated how 
the equation set may be reduced to scalar equations for the case of scalar perturba- 
tions in a K = universe. The extension to K ^ and vector or tensor modes is 
straightforward. The results of a numerical solution of the equations for scalar per- 
turbations will be given elsewhere 0. We expect that these calculations will confirm 
the result of other groups who have made their calculations by working carefully in 
specific gauges. 

4 CMB Anisotropics in Global Defect Theories 

Our final topic is a review of the results of some recent calculations by Pen, Seljak 
and Turok || [7j of the power spectra in global defect theories. Unlike the calculation 
of CMB anisotropies in inflationary models, the calculation for defect theories has 
been plagued by computational difficulties. These difficulties arise from the contin- 
ual generation of perturbations by the causal source (the defects). One must solve 
for the non-linear evolution of the source, as well as the linearised response of the 
Einstein/fluid/Boltzmann equations to the causal sources. 

Recently, an efficient technique has been developed for the calculation of power 
spectra in defect theories ||, which makes use of the unequal time correlator of the 
defect stress-energy tensor. This technique has provided the first accurate computa- 
tions of the CMB and matter power spectra in global defect theories. The results of 
the calculation for the CMB power spectra (the Cj) are shown in Figure [|, for global 
defects comprised of strings, monopoles, textures and non-topological textures. The 
contributions from scalar, vector and tensor modes are shown separately, as well as 
their total. The most striking feature of the figure is the large contribution from 
vector modes to the power on large angular scales. The vector modes dominate the 
total signal for / < 100, at which point they are suppressed by the horizon size at re- 
combination. Vector modes are insignificant in inflationary models since the vorticity 
decays as the universe expands to conserve angular momentum. However, in defect 
theories, the defects are a continual source of vorticity giving rise to a significant 
vector component to the CMB power spectrum. In addition, the continual (causal) 
sourcing of the fluid perturbations leads to decoherence, where scalar modes of a given 
scale (inside the sound horizon) oscillate with different phases in different parts of the 
universe, compared to the coherent oscillations which occur in inflationary models. 
This decoherence manifests itself as the suppression of the secondary Doppler peaks 
in Figure |j. 
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Figure 4: The contributions to the CMB power spectrum from scalar, vector and 
tensor components in global defect models with strings, monopoles, textures and 
non-topological textures. Reproduced with permission from Pen, Seljak and Turok 
(1997). 

The significant contribution of the vector modes to the large scale power in the 
CMB spectrum leads to a suppression in the total signal for I > 100 compared to 
inflationary models, when both are normalised at large angular scales. In Figure [5] 
we show a comparison of the predictions of global defect theories, with COBE nor- 
malisation at I = 10, to the current generation of CMB measurements. The effect of 
vector modes and decoherence is to leave the predictions of the defect models system- 
atically below the current degree-scale data. A similar effect occurs with the matter 
power spectrum, P(k). The matter power spectrum calculated in global defect theo- 
ries is compared with the predictions of a standard cold dark matter model (sCDM, 
with Hq = 50kms _1 Mpc -1 and a scale-invariant spectrum of initial perturbations) 
in Figure |6|. The defect calculations and the sCDM calculations are normalised to 
COBE, and are compared with the matter power spectrum inferred from the galaxy 
distribution. The large contribution of vector modes to the CMB power spectrum in 
the range of the COBE normalisation, suppresses the scalar component compared to 
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Figure 5: A comparison of the CMB power spectra calculated in defect models with 
the current observations. The calculated curves are normalised to COBE at / = 10. 
Reproduced with permission from Pen, Seljak and Turok (1997). 



sCDM. The COBE normalised defect models predict too little power in the matter 
spectrum on larger scales, and too much power on smaller scales. As is well known, 
the sCDM model fares much better on large scales, but also predicts too much small 
scale power. 

To conclude, current degree-scale observations of the CMB power spectrum, and 
the matter power spectrum inferred from galaxy clustering, do not favour global 
defect theories of structure formation, on the basis of the first accurate computations 
of the spectra in these theories. The significant contribution to large scale power from 
vector modes, which occurs in all defect models, leads to systematically low (COBE 
normalised) predictions for degree-scale CMB observations, and for the matter power 
spectrum. 
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Figure 6: The matter power spectra computed in COBE normalised defect theories 
(lower curves) compared with the predictions of the standard CDM (sCDM) model 
(upper curve) for the same normalisation. The data points show the power spectrum 
inferred from galaxy distributions. Reproduced with permission from Pen, Seljak and 
Turok (1997). 
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